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ABSTRACT 


This paper investigates a class of decomposable codes, their distance and structural 
properties. It is shown that this class includes several classes of well known and efficient 
codes as subclasses. Several methods for constructing decomposable codes or decomposing 
codes are presented. A two-stage soft decision decoding scheme for decomposable codes, 
their translates or unions of translates is devised. This two-stage soft-decision decoding is 
suboptimum, and provides an excellent trade-off between the error performance and decoding 
complexity for codes of moderate and long block length. 
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1. Introduction 

To decode a long block code with maximum likelihood decoding is practically impossible 
because the decoder complexity is simply enormous, However, if a code can be decom- 
posed into constituent codes with smaller dimension and simpler structure, it 
is possible to devise a practical and yet efficient scheme to decode the code as follows: (1) 
The constituent codes are decoded sequentially in multiple stages. The decoded codeword 
of one constituent code at one stage is passed to the next stage for decoding the next con- 
stituent code; and (2) From the decoded constituent codewords, we form the overall decoded 
codeword based on the code structure. This type of decoding is not optimum even if the 
decoding of each constituent code is optimum. It is a suboptimum decoding scheme. How- 
ever, this decoding scheme reduces the overall decoding complexity drastically comparing 
with the optimum decoding of the overall code. This is because the constituent codes are 
smaller in dimension and much simpler in decoding complexity. If the constituent codes 
have the right structure (such as trellis structure ) and the decoding scheme is devised 
properly, an excellent trade-off between the error performance and decoding complexity 
can be attained. 

In this paper, we first define a class of decomposable codes, and show that several classes 
of good known codes are decomposable. Several methods for constructing decomposable or 
decomposing codes are presented. Then a two-stage suboptimum decoding scheme is devised 
for decomposable codes, their translates and unions of their translates. This suboptimum 
decoding scheme reduces the overall decoding complexity drastically comparing with the 
single-stage optimum decoding for the same decomposable code while maintains excellent 
performance. The error performance of some specific decomposable codes based on the 
proposed two-stage suboptimum decoding is evaluated and simulated. It is shown that the 
proposed two-stage suboptimum decoding achieves excellent error performance with reduced 
decoding complexity. 

2. Decomposable Codes 

Let I be a finite set of symbols on which an addition “+”, a substraction 
and a distance measure d(-, ■) (Euclidean or Hamming) axe defined. Let 0 denote the 
zero element of L . The distance measure d(u,v) between two elements, u and v, in L 
is assumed to satisfy the following properties: (1) d(tz,u) = 0 if and only if u = v, and 
(2) d(u,v) — d(v,u). For the binary case, L is defined as the binary set {0, 1}, the distance 

measure d is given by d( 0, 1) = 1, and the on L is defined as the modulo-2 addition. For 
a 2*- ary PSK signal set 5, L is simply defined as the set of all binary label strings of length 
i for the signal points in 5. In this case, the addition u +” means the bit-wise modulo-2 
addition. 

For two ./-tuples, u = (u 1? u 2 , • • ■ ,u„) and v = (t>i, t> 2 , • • • , v„) over L, let u + v and d( u, v) 
be defined as, 

u + v = (til + *h,iz 2 + t>2, • • • ,Uj + Vj) (2.1) 
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4“.v) = 52 </(«,, U.) ( 2 - 2 ) 

»=1 

respectively. For a block code C over L, let d[C] denote the minimum distance of C with 
respect to the distance measure d. 

For block codes C and C' of the same length over L, let C + C‘ denote the set {u + v : 
u € Cand v 6 C'}. If C' consists of a single codeword v, then C + {v} is called a translate 
of C [1]. If for any two different codewords u and u' in C and any two different codeword 
v and v ' 6 C', u + v ^ u' 4- v',then C and C' axe said to be “independent”. If a code 
is closed under the component-wise addition “+” and subtraction the code is said to 

be additive. For an additive code C and its supercode C' which is a set of cosets of C, 
let [ C'/C ] denote the set of representatives of cosets of C in C'. For C itself, the all-zero 
n-tuple 0 is always chosen as its representative. 

Let m be a positive integer, and let n be a positive inetger divisible by m. For a code C 
of length n over L, C is said to be decomposable with respect to U x and U 2 if there exist 
independent codes U x and U 2 of length m over L, and codes C\ and C 2 of length n over L 
such that : (Al) C x C (U x ) n ' m (A2) C 2 C {U 2 ) n ' m , and (A3) C = C x + C 2 = {u + v : u € 
Ci and v € C 2 }, where A"/"' denotes the following Cartesian product of A, 

A n / m = A x A x ■ • • x A . 

' v " 

n/m 


For u and u' in U x and v and v' in U 2 , let dfj/uju, u') and dff ^Jv, v') be defined as follows: 

i u') = min{d(u + v, u' + v') : v and v' in U 2 } (2.3) 

) = min {^( u + v,u + v') : u in U x } (2.4) 

Clearly ^^(u, u) = d^(v,v) = 0, ^^(u, u') = ^^(u', u) and v ') = 

dvi,u 2 ( v » v )* 

For simpicity, we regard an n-tuple over L as a n/m tuple over L m , the set of all m-tuple 
over L, and vice versa. We also consider a code of length n over L as a code of length n / m over 
L m and vice versa. For two codewords , x = (xi,x 2 , • ■ • ,x„/ m ) and y = (yi,yj, • • • ,y n /m) i n 
Ci with x } and y } in U, for 1 < i < 2 and 1 < j < n/m , define the distance between x and 

y based on the distance measure ,{/,(-, ■) given by (2.3) or (2.4) as follows: 

n/m 

d { ul lh (x,y) = ( 2 -5) 

j = i 

Then the minimum distance of C, based on the distance measure given by (2.3) or (2.4) is 
defined as follows: 

dv t ,u 2 [ c i\ = {^S,t/ 2 ( x >y) : x >y € C,}. (2.6) 
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Lemma 1: The minimum distance D[C ] of a decomposable code C(= C x -f C 2 ) with respect 
to U\ and U 2 is lower bounded by : 

i[C] > 4?,VJCj]). (2.7) 

If the equality in (2.7) holds, then C is said to be strictly decomposable. 

It is clear that given two block codes C x and C 2 over L which satisfy the conditions (Al) 
and (A2), we can form a decomposable code C = C\ + C 2 over L. In the following, we use 
some examples to demonstrate how to construct decomposable codes and to show that some 
known codes are decomposable. Hereafter, the following notations are used: (1) V„ denotes 
the set of all binary n-tuples (= {0,1}" ) ; (2) P„ denotes the (n,n - 1) linear code which 
consists of all the even- weight binary n-tuples; (3) RM, :} denotes the j-th order Reed-Muller 
(RM) code of length n = 2' ; and (4) For a code C, C x denotes the dual code of C. P x 
consists of the all- zero n- tuple 0 and the all-one n-tuple 1. 

Example 1 : Consider a binary case where m = 2 , U\ = {(0, 0), (0, 1)}, U 2 = P% — 
{(0,0), (1,1)} and U l + U 2 = {0,1} 2 . Then <,^((0,0), (0, 1)) = <*((0,0), (0, 1)) = d(0,l) 
and ^,((0,0), (1,1)) = °)» (*» *)) = 2<*(0,1). Let n be an even positive integer. 

For a binary nj 2- tuple, u = (u l5 u 2 , ••• u n/2 ), let <£i(u) and <£ 2 (u) denote two n-tuples 
obtained from u as follows : (1) <£i(u) is obtained from u by substituting 00 and 01 for 
0 and 1 in u respectively. (2) ^ 2 (u) is obtained from u by substituting 00 and 11 for 0 
and 1 in u respectively. Let X and Y be two linear binary block codes of length nj 2. 

Define the following two linear binary codes of length n : C\ = <^i[X] = {<f> i(x) : x E X} 
and C 2 = <j> 2 [Y] = {^ 2 ( y) : y G Y}. Then C = C x + C 2 is decomposable with respect to 
Ui = {(0,0), (0, 1)} and U 2 = {(0,0), (1, 1)}. The above construction of a decomposable code 
with respect to U x = {(0,0), (0, 1)} and U 2 = {(0, 0), (1, 1)} is simply a permutation of the 
|u|u + v|-construction of codes. Therefore, codes obtained from the |u|u + v|-construction 
are decomposable with respect to Ui = {(0,0), (0, 1)} and U 2 = {(0,0), (1, 1)} after a certain 
permutation. 

Now we show that RM codes fire decomposable with respect to = {(0,0), (0, 1)} and 
U 2 = {(0,0), (1,1)}. Let Ci = <f>i[RM m -i t r-i] and C 2 = <j> 2 [RM m - lir }- It is known [1] that 

RMms = <j>\[RM m -\ t r-\] + URM m -lA- (2-8) 

From (2.3) and (2.4), we find that d[^ f r 2 [Ci] = 2 m r and d[/ l ^ l 7 J [C 2 ] = 2 m T . Since 
d(RM mjr ) = 2 m ~ r , the equality of (2.7) holds. Therefore Reed-Muller codes are decom- 
posable with respect to U] = {(0, 0), (0, 1)} and U 2 = {(0,0), (1, 1)}. 

Example 2 : Consider a binary case where m = 4, U n = {(0, 0), (0, l)} 2 , U i2 = 

{(0,0, 0,0), (0,0, 1,1)}, U x = U n + U u and U 2 = P 4 X . Let C' n be a block binary code 
of length n/2 and and minimum Hamming distance C' u be a binary block code of length 
n/4 and minimum Hamming distance <5 i 2 , and C' 2 be a binary block code of length n/4 
and minimum Hamming distance S 2 . For a binary j-tuple u, let ^ 2 <^i(u) ( or ^ 2 <^ 2 (u) ) 
denote the 4.;-tuple derived from u by substituting (0,0, 0,0) for each component 0 and 
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(0,0, 1, 1) (or (1, 1,1,1)) for each component 1. Let C n ,C n ,C x and C 2 be defined as follows 
: (1) C n = <t>i[C' u }\ (2) C n = M i[C[ 2 ]\ (3) Cj = C„ + C 12 ; and (4) C 2 = <hfa[C' 2 \, where 

<j> x and </> 2 are defined in Example 1, and for a mapping / and a block code C , f[C ] = 
{/(v) : v € C'}. Then C = C, -f C 2 is decomposable with respect to JJ\ and I/ 2 . It can be 
shown that d[C] > mm{6 n ,26i 2 ,4£ 2 }. If C‘ U ,C' U and C\ contain the all-zero n-tuple 0, then 
d[C] = min{<5 n ,2<5 12 ,4<5 2 }, and C is strictly decomposable. 

As an example, consider RM m<r with r > 2. From (2.8), we see that RM m , r is decom- 
posable with respect to U\ and U 2 as follows: 

RM m>r = {<f>x[RM m - l,r-l] + 

-f-^> 2 (^ 2 [i?A/ m _ 2l r]. (2-9) 

In fact, C\ is decomposable with respect to U\\ + U\ 2 . 

Besides RM codes, multi-level concatenated codes (2] and multilevel modulation codes 
[3-5] are also decomposable. Some primitive BCH codes are union of a decomposable code 
and its translates. 

3. Two-Stage Decoding 

In this section, we present a two-stage suboptimum decoding for a decomposable code, 
Q = C\ + C* 2 , over L with respect to U\ and U 2 which achieves the distance min{d[^ u 7 \C\) 1 

du! u 2 [^ 2 ] } • The first stage is for decoding C x based on the distance measure d[^ {/j (-, •), and 
the second stage decoding is for decoding C 2 based on the distance measure <4/^1/, (*, •)• 

Let C = C\ -f C 2 be a decomposable code of length n over L with respect to U\ and 
U 2 over L m , where n is divisible by m. Suppose C is used for error control over an AWGN 
channel. Assume that all the codewords in C are equally likely to be transmitted. For 
v E L, let s(u) denote the signal point in R h represented by v, where R h denotes the set 
of all ^-tuples of real numbers, and for a j-tuple v = (iq,u 2 , • ■ • ,Vj) over L , let s(v) denote 
the tuple, (s(ui), s(u 2 ), • • • , ^(uj)). For 2 and z in R h , let || 2 - z || denote the Euclidean 
distance between z and z . For ^-tuples z = (rj, 2 2 , • • • , Zj) and z = (zj, z 2 , • ■ ■ , Zj) over R h , 
let || z - z' || 2 be defined as Ei=i II z < “ z \ II the squared Euclidean distance between z 
and z\ The distance measure d on L is assumed to satisfy the condition that for u and v 
in L, || s(u) - s(u) || 2 > d(u,v). Suppose that (1) for u € L, a one-to-one mapping T u from 
R h to R h which preserves the squared Euclidean distance is defined, and (2) for u and v in 
L, s(v — u) = T u (s(i’)). Suppose a codeword in C is transmitted. Let z = (zi, z 2 , • ■ • , ^n) f> e 
the received vector over R h . A two-stage decoding procedure (D) for C can be formulated 
as follows: 

(Dl) Decode z into a codeword in C\. Let u = (u x , u 2 , ■ ■ • , u n ) be the decoded codeword 
in Ci. 

(D2) From u, we form the following vector based on the mapping T u (-): 

T u (z) = (T Ul (=i),T U3 (z 2 ),---,T Un {z n )) 
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Decode !T u (z) into a codeword in C 2 . Let v= (t>j, v 2 , • • ■ > v„) be the decoded 
codeword . 

(D3) The decoded codeword is given by u + v = («i + t>i, u 2 + u 2 , • ■ • » *<« + v„). 


The decodings of Ci and C 2 may be either soft-decision or hard-decision, maximum 
likelihood or bounded distance. We may take advantage of the structure of Cj, C 2 , U\ 
and U 2 to simplify the decodings. Since the dimensions of C\ and C 2 are smaller than that 
of C and U x and U 2 are much shorter than C, the above two-step decoding may result in a 
drastic reduction in decoding complexity . 

Suppose both C\ and C 2 have ?i/?n-section trellis diagrams, each branch in the trellis 
diagram for is a codeword in U\ and each branch m the trellis diagram for C 2 is a 
codeword in U 2 . Then a soft-decision (SD) decoding procedure for the above two-stage 
decoding can be devised as follows: 


(SDl) Divide the received vector z into n/m sections, z ,, for 1 < i < n/m, with each 
section containing m consecutive components of z. Using the trellis diagram for 
C\ and the Viterbi algorithm, decode the n/m-tuple z = (zj, z 2 , • • • , z n / m ) into a 
codeword u in C\ . The branch metrics in the i-th section of the trellis diagram 


are 



(3-1) 


for each r € U\. If U 2 has a trellis diagram, the branch metric of (3.1) can be 
computed by using a trellis diagram for {r} + U 2 . 


(SD2) Let u = (ui,u 2 , ■■■,«„) be the decoded codeword in Ci at step SDl . Form the 
vector 

T u (z) = (T u , ( 21 ), T U j(2 2 )j ' ■ ' j T Un ( 2 n )), 


and decode it into a codeword u in C 2 based on the trellis diagram for C 2 using 
the Viterbi algorithm. To compute the branch metrics at each trellis section, we 
divide T u (z) into n/m sections, T u (z,) for 1 < i < n/m. The branch metrics in the 
i-th section of the trellis diagram for C 2 are computed based on ||T u (zj) - s(v)|| 
with v € U 2 . 


(SD3) The decoded codeword in C is given by u + v. 


It is clear that if C 2 is decomposable, C 2 can be decoded in two stages. The decoding 
procedure (D) can be applied recursively. 

Note that the two-stage soft-decision maximum likelihood decoding is a suboptimum de- 
coding. In the next section, we will show that this suboptimum decoding achieves excellent 
error performance comparing with the single-stage optimum (soft-decision maximum likeli- 
hood ) decoding while reduces the overall decoding complexity drastically. It is clear that a 
two-stage hard-decision decoding can be devised. 
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The proposed two-stage decoding is applicable to any code which is a translate of a 
decomposable code. If a code C is a set of translates of a decomposable code, the decoding 
can be carried out as follows: (1) Decode the received vector z into a codeword in each 
translate by the above two-stage decoding ; and (2) Choose the most probable one among 
all the decoded words as the final decoded codeword. 

Sufficient conditions under which the proposed two-stage decoding procedure result in a 
correct decoding for an AWGN channel, are given in Lemma 2. 

Lemma 2 : Using the two-stage soft decision maximum likelihood decoding, if a codeword 

w € C is transmitted and z is received such that 

j| z — s(w) || 2 < min{d[^c; J [Ci] , [C 2 ] } /4 

then z is correctly decoded into w. 


4. Performance 

In the following, we use three examples to demonstrate that the proposed decoding of 
decomposable codes indeed reduces the decoding complexity drastically while maintains 
excellent error performance. A decomposable code may be decomposed in different ways. 
Different decompositions result in different error performance and decoding complexities. In 
general, the error performance of the first stage decoding has dominant effect on the overall 
error performance and the first constituent code should be chosen as large as possible within 
allowable decoding complexity. 


Example 3. : In Example 1, we showed that the r-th order RM code C = RM miT 

of length 2 m and minimum distance 2 m-r is strictly decomposable with respect to U x = 
{(0,0), (0, 1)} and U 2 = {(0, 0), (1, 1)}. The constituent codes of C = RM m , r are: Cj = 
r -i] and C 2 = <j> 2 [RM m -\,r]- RM m - i.r-i and RM m . hr have 4-section trellis 

diagrams with 2(™-*) states and 2 C”'*) states respectively^]. It is easy to see that C x and 
C 2 also have 4-section trellis diagrams with 2( r - 1 ) states and 2^ ■■ ) states respectively. The 
RM code C = RM m<r has a 4-section trellis with 2^ >■ *) states. If we decode C = RM m ,r with 
the optimum decoding using Viterbi algorithm, the state complexity of the Viterbi decoder 


is 2^ m r *). However, if we decode C = RM m<r with the proposed two-stage suboptimum 
decoding algorithm, the overall state complexity of the Viterbi decoder is 2^ «• ) + 2^ r - 1 ) 
which is much smaller than 2( " r \ the state complexity of the 4-section trellis diagram for 


the overall code C = for m > 5. For example, let vn — 6 and r — 2. Then C — 

is a (64,22) code which has a 4-section trellis diagram with 2(‘) = 1024 states. However, its 
constituent codes, C x - and C 2 = <j> 2 [RM$, 2 ) , have 4-section trellis diagrams with 


2 ( 1 ) = 16 states and 2^) = 64 states respectively. Using the two-stage suboptimum decoding, 


the total state complexity of a Viterbi decoder is 80 which is much smaller than the state 
complexity 1024 for the single-stage optimum decoding for C = RM 6<2 . The code C = RM 6<2 
has minimum distance 16, which is comparable to a rate - 1/3 optimum convolutional code 


of constraint length 8 and free distance 16 [7]. To decode this convolutional code with the 
Viterbi algorithm, a Viterbi decoder of 128 states is needed. 
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The error performance of the RM code C = RM &, 2 with the proposed two-stage sub- 
optimum soft-decision decoding based on the above decomposition, denoted 2si, is shown 
in Figure 1. We see that it achieves 5.1db coding gain over the uncoded BPSK at block 
error rate 10 -5 . Figure 1 also shows simulation results pj^ f [2si] on the probability that the 
first stage decoding is correct but the second stage decoding is incorrect. Compared with 
the simulation results p ic , 4 [2si] on the overall error probability, is very small. This 

fact suggests to choose a supercode of <f> 1 [ RM54] as the first constituent if the decoding 
complexity is allowable. By applying the decomposition shown in Example 2 to RM 6l 2, we 

have Cj = 4>i[RMs,i] + <t>-[4>2[RkU,i] and C 2 = fafalRhUa] as the first and second constituent 
codes, respectively. C\ and C 2 have 4-section trellis diagrams with 128 states and 8 states, 
respectively. Note that the first consituent code in the decomposition contains the first con- 
stituent code of the former decomposition as a small subcode. The error performance of 
RMq <2 with this decomposition, denoted 2s 2 ) is shown in Figure 1. It achieves 5.6dB coding 
gain over the uncoded BPSK at block error rate 10 -5 . Figure 1 also shows simulation results 
, denoted p iCi ,[l-s] , and union upper bounds, denoted p ^ [Is], on the probability of incor- 
rect decoding for the single stage optimum decoding. We see that the coding losses of the 
two-stage suboptimum decodings based on decompositions 2sj and 2s 2 from the single-stage 
optimum decoding are about l.OdB and 0.5dB, respectively, at the block error rate 10 -5 . 
However there is a big reduction in decoding complexity. 

It is known that an extended and permuted BCH code of length 2 m and minimum distance 
2 m-r contains the r-th order Reed-Muller code RM m>r as a proper subcode[l,7]. As a result, 
this extended and permuted BCH code is a union of cosets (or translates) of the Reed-Muller 
code RM m , r . Since RM m , r is decomposable and has relatively simple trellis structure, this 
extended and permuted BCH code can be decoded with the proposed suboptimum soft- 
decision two-stage decoding. First we decode the received n-tuple z into a codeword in each 
translate of the RM m , r code by the proposed two-stage decoding. Then we choose the most 
probable one among all the decoded words as the final decoded codeword. All the translates 
of the RM m ,r code have trellis diagrams isomorphic to that of the RM m<r code. Let K be 
the total number of translates of RM m>r ( including RM m<r itself). If A is not too big, it is 
possible to implement K separate but identical Viterbi decoders to decode the K translates 
of RM mt r in parallel. This definitely speeds up the decoding process. 

Example 4: Let C be the extended and permuted (64, 24) BCH code whose complexity 
of trellis structure has been analysed in [8]. C contains the (64, 22) second order Reed-Muller 
code as a proper subcode. There are 4 cosets in C modulo RM§ t2 . The coset code \C j RM s i2 ] 
is generated by the two codewords gi and g2 which can be determined [8]. Each coset of 
C modulo RM 6i2 can be decomposed into two constituent codes, C\ and C2, either in the 
form given in Example 1 ( denoted 2sj ) or in the form given in Example 2 ( denoted 2s 2 ). 
The two constituent codes, Ci and C 2 , have a 4-section trellis diagrams with 16 states and 
64 states or 128 states and 8 states, respectively. Hence each coset can be decoded with the 
suboptimum soft-decision two-stage decoding with two Viterbi decoders, one with 16 states 
(or 128 states ) and the other with 64 states ( or 8 states ). The overall decoder for the 
(64,24) extended and permuted BCH code consists of 4 separate two-stage decoders. These 
four decoders process the received 64-tuple in parallel. The total number of states for the 
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overall decoder is 320 (or 544). The error performance of this BCH code with suboptimum 
soft-decision two-stage decoding based on the decomposition 2si is simulated and shown in 
Figure 2. We see that it achieves 5.4 dB coding gain over the uncoded BPSK at block 
error rate 10 -5 , which compares favorably with the rate-1/3 Odenwalter convolutional code 
of constraint length 8 and free distance 16. If each coset of the (64,24) extended BCH code 
modulo RMe ,2 is decoded with one-stage optimum decoding, four 1024-state Viterbi decoders 
are needed, one for each coset . This amounts to a total of 4096 states which is much greater 
than the number of states ( 320 or 544 ) in a two stage decoding as described above. 

Example 5 : Let C be the extended and permuted (64, 45) BCH code whose complexity 
of trellis structure has been analyzed in [8]. C includes RMe, 3 , a (64,42) code, as a subcode. 
[C/RMej) consists of eight cosets. Each coset can be decomposed into two constituent codes, 
Ci and C 2 have 4-section trellis diagrams with 64 states and 16 states, respectively. Hence 
each coset can be decoded with the suboptimum soft-decision two-stage decoding with two 
Viterbi decoders, one with 16 states and the other with 64 states. The overall decoder for 
the (64,45) extended and permuted BCH code consists of 8 separate two-stage decoders. 
These 8 decoders process the received 64- tuple in parallel. The total number of states for 
the overall decoder is 640. The error performance of the BCH code with suboptimum soft- 
decision two-stage decoding based on the above decomposition, denoted 2sj, is simulated 
and shown in Figure 3. We see that it achieves 5.4dB coding gain over the uncoded BPSK 
at block error rate 10 -5 . If each coset of the (64,45) extended BCH code modulo RM^z 
is decoded with one-stage optimum decoding, eight 1024-state Viterbi decoders are needed, 
one for each coset. In this case, the total number of states is 8192 which is much greater 
than that of the two-stage decoding described above ( 640 states ). 

From the above three examples, we see that the proposed code decomposition and two- 
stage suboptimum decoding achieves excellent reduction in decoding complexity. 

5. Conclusion 

In this paper, we have shown that code decomposition with multi-stage suboptimum de- 
coding achieves excellent error performance with a drastic reduction in decoding complexity. 
This finding may have an impact on the future designs of error control systems for reliable 
data communications. For future work, we should focus in searching efficient decomposable 
codes with simple trellis structure so that soft-decision Viterbi decoding can be applied. 
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Figure 1. Error performance of RM 6>2 
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